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We introduce a self-consistent stochastic coarse-graining method, which includes both metric 
and scalar field fluctuations, to investigate the back reaction of long wavelength perturbations in 
single-scalar driven inflation, up to the second (one loop) order. We demonstrate that, although back 
reaction cannot be significant during the last 70 e-foldings of inflation with a smooth potential, there 
exist non-smooth inflaton potentials which allow significant back reaction, and are also consistent 
with cosmological observations. Such non-smooth potentials may lead to the generation of massive 
primordial black holes, which could be further used to constrain/verify these models. 



PACS numbers: 98.80.Cq, 98.80.Qc 



I. INTRODUCTION 



In the context of inflationary cosmology, initial seeds 
of today's structures are generated in the de-Sitter phase 
of an inflationary universe, as a consequence of quantum 
vacuum fluctuations (see e.g. p| for a comprehensive 
overview of the theory of cosmological fluctuations and 
for a recent introductory overview). The standard 
way of analyzing these fluctuations is through linear per- 
turbation theory. However, the linear analyses is limited 
by the non-linear nature of the Einstein equations, and in 
particular, the presence of perturbations is likely to affect 
the evolution of the background cosmology at the non- 
linear level. The back reaction of short wavelength grav- 
ity waves on an expanding Friedmann- Robertson- Walker 
(FRW) cosmology is a well- understood problem Q ■ How- 
ever, in the context of inflationary and post-inflationary 
cosmology, the scalar metric fluctuations (fluctuations 
coupled to energy density and pressure perturbations) 
are believed to dominate over the effects of gravity waves, 
and thus, it is of great interest to understand the possible 
back reaction of these scalar perturbations on the cosmo- 
logical background. Furthermore, during inflation, the 
phase space of infrared modes (defined as modes with 
wavelength greater than the Hubble radius) grows ex- 
ponentially, whereas the phase space of the ultraviolet 
modes does not grow and since the amplitude of the as- 
sociated metric fluctuations of these infrared modes does 
not decrease in time , the back reaction of these infrared 
modes may grow to be significant. There have been dif- 
ferent approaches to address this question. In 0, @, 
the effective energy-momentum tensor formalism of [3| 
was generalized to study the back reaction of infrared 
modes on the spatially homogeneous component of the 
metric. The result was that, in a slow-roll inflationary 
background, the back reaction of infrared modes takes 
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the form of a negative cosmological constant. This was 
later confirmed in Q using very different techniques. 

The objections (see 0) to these analyses are mainly 
concerned with the fact that 0, l| do not consider a 
local observable quantity, and even if the effect is not 
simply a gauge artifact, it is still not clear how a lo- 
cal observer could distinguish the back reaction of long- 
wavelength modes from an initial condition ambiguity in 
a homogenous universe. More importantly, calculating 
an "observable" from the spatially averaged metric will 
not in general give the same result as calculating the spa- 
tially averaged value of the observable. A later work by 
Abramo and Woodard in Q tried to resolve this problem 
by identifying a local physical variable which describes 
expansion rate of the universe, then calculating the back 
reaction of cosmological perturbations on this quantity. 
A different approach was introduced in , where a sim- 
ple variable describing the local expansion rate was con- 
sidered and calculated up to second order in the metric 
fluctuations in a model with a single inflation field [28| . In 
this approach, the leading contribution of back reaction 
of infrared modes to this variable was evaluated. When 
evaluated at a fixed value of the matter field ip, the only 
physical clock available in this simple system, the depen- 
dence of the expansion rate on the clock time takes on 
exactly the same form as in an unperturbed background 
pfjf . Thus, the leading infrared back reaction terms had 
no locally measurable effect in this system. Similar con- 
clusions were reached in 

However, the fact that during inflation the phase space 
of infrared modes is growing in time, i.e. modes are ex- 
iting the Hubble radius, leads us to the conjecture that 
there might be a locally measurable effect just due to 
these crossings, and if we include the corrections due 
this effect, the results for single scalar field models may 
change. Since these modes are generated in the quan- 
tum mechanical ultraviolet limit, handling them in gen- 
eral is very difficult. Attempts, so far, in order to model 
the evolution through writing a classical Langevin-like 
equation for the coarse-grained scaler field, the so-called 
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stochastic approach |l3[ . fail to include the inflaton and 
metric fluctuations self-consistently |3(|. In this work, 
for the first time, we devise a self-consistent stochastic 
framework to study the evolution of the coarse-grained 
inflaton. While the approach uses the results of linear 
(metric+inflaton) perturbation theory for the ultravio- 
let and transitional modes, it is rigorously extended to 
the non-linear infrared regime, which allows a consistent 
study of the evolution of fluctuations through their entire 
history. We then apply this framework to the problem 
of back reaction in inflationary cosmology, and investi- 
gate how/when the non- linear effects may change prop- 
erties/predictions of the standard slow- roll inflation. 

The structure of the paper is as follows: In Section 
II, we explain how the linear theory of perturbations can 
be self-consistently generalized to the non-linear infrared 
regime, including both inflaton and metric fluctuations. 
We also introduce 0, the local expansion rate, and how 
it is related to metric perturbations. In Sec. Ill, we show 
that the quantum generation of fluctuations can be mod- 
elled as a stochastic term in the equation of motion for 
O, coarse-grained over local Hubble patches. In Sec. IV, 
we apply this equation to a chaotic inflationary scenario 
with a quadratic potential V(ip) = \m 2 ip 2 and evaluate 
the local expansion rate with back reaction corrections 
due to infrared modes. Sec. V contains the generaliza- 
tion of the result from Sec. IV to any power-law poten- 
tial. In Sec. VI, we investigate if back reaction can ever 
significantly change the inflationary predictions for our 
observable universe. Finally, Sec. VII summarizes our 
results and concludes the paper. Throughout this paper, 
we use the natural units, and set the Planck mass equal 
to unity, i.e. 8ttG = c = h = 1. 



II. GENERALIZED BARDEEN PARAMETER 
IN THE NONLINEAR REGIME 

Our analysis is based on the consistent quantization of 
the linearized metric and matter perturbations in a clas- 
sical expanding background space-time. In this frame- 
work, it can be shown that, in the linear regime, the 
Sasaki-Mukhanov parameter v 0, , defined as yl| : 



We use the following ansatz for our non-linear metric: 



ds 2 



z 2<t >dt 2 



(4) 



which is a generalization of a metric with linear scalar 
perturbations. This ansatz, though not a general solu- 
tion, represents the most important sector of the met- 
ric, as any such metric can be continuously connected to 
a metric with linear scalar perturbations, which is the 
dominant sector if perturbations are generated during 
slow-roll inflation [Tj|. We consider the simplest model 
of matter with a single scalar field ip, the inflaton. We 
also use a perfect fluid model to describe the expansion 
of the universe. The velocity four-vector field u a can be 
considered to be tangential to a family of world lines of 
a matter fluid in a general space-time. This four-vector 
is normalized such that 



u a u a = 1 , 



(5) 



where a runs over the space-time indices. In terms of 
this four-vector, the local expansion rate is defined by 



= 



= 3ff, 



local i 



(6) 



where H\ oca \ is the local Hubble constant. It represents 
the local expansion rate of the tangential surfaces orthog- 
onal to the fluid flow. It is easy to verify that, for our 
specific choice of gauge, follows the following relation 







(7) 



where a prime denotes the derivative with respect to r, 
the proper time, and is related to the coordinate time via 
dr = e^dt. The Lagrangian and the energy-momentum 
tensor for the inflaton are given by: 



c = ±d"<pd a <p-v{<p), 
T2 = d»<pdu<p - est;. 



(8) 
(9) 



It can be shown that for slow-roll inflation, in the lone 
wavelength limit, satisfies the following condition |9( 



v = a[S(p + ( 



H 



axp 

IT 



C, 



(i) 



obeys the equation of motion of a free scalar field, on 
scales much smaller than the Hubble radius. Here, £ is 
the Bardeen parameter |l7j : 



C - cj>-^-{H<t> + 4>), 

ti 



(2) 



where <p is the scalar metric perturbations in the longi- 
tudinal gauge, in a spatially flat universe: 



ds 1 



(1 + 2(j))dt 2 - a 2 (t)(l - 2(/))Si j dx i dx j 



(3) 



= y/3V(<p). 



(10) 



In [19| it is demonstrated that, in these limits, the gen- 
eralized Bardeen parameter: 



c(x) = V ' 2 J { ^ P ~ ) d ^ 



(11) 



remains constant in the classical non-linear regime |20j . 
Here, g is a function of the inflaton field, obtained by 
integrating the Goi Einstein constraint equations, and 
satisfies 



2iP' = g(fp); ^ 



dip 



(12) 
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It is easy to see that for linear perturbations S((x) re- 
duces to the usual Bardeen parameter £, which is known 
to remain constant for adiabatic super-Hubble modes 

Combining the above equation, with Eqs. Q and ( 
II 111 . we see that 



~, e 1 aiog(e) x , 



(13) 



in the classical long- wavelength limit. In the next section, 
we show how quantum fluctuation affect Eq. i|13|) . 



III. COARSE-GRAINING THE EXPANSION 
RATE 

In this section, we use coarse-graining on the scale of 
the Hubble radius to provide a framework to analyze back 
reaction, and obtain a new derivation of the equation 
of motion for stochastic inflation [2l|, which takes into 
account the metric fluctuations. 

Let us define the coarse-grained generalized Bardeen 
parameter, £ c , in the Fourier space, as: 



Cc(k) = 0(aH - fc)C(k) 



(14) 



where 9 is the step function (not to be mistaken with 
©, the local expansion rate) which is zero for k > aH, 
and is one for k < aH |34| . Taking the time derivative 
of the above equation; the time derivative of £ vanishes 
for infrared modes and the time derivative of 9 is a delta 
function at k = aH, which yields 



£(k) = c(k)(fr + |oa(i-^). (i5) 

This simply implies that C' c ( x ) m rea, l space is related to 
C{k) in the following way 

e ( x, = ( ff+ |)/^e<'Wl-A ) . (16) 

Two-point function of £(fc) can be described in terms of 
its power spectrum, P^(k), defined by: 

(C(ki)C(ka)) = (2^) 3 ,5(k 1 - k 2 )P f (fe). (17) 

Assuming that Hubble-scale perturbations are linear, 
P^{k) can be derived in the linear theory of quantum 
generation of perturbations , 



H 4 



(18) 



Using this knowledge for £(&) and assuming H' <C H 2 
(the slow-roll condition) in the above equation, we can 
calculate the two-point function for (' c (x, r). Note that 



H in Eqs. I|14I18|I is the local Hubble constant, and can 
be replaced by 9/3. Combining Eqs. I|16ll8f) . after per- 
forming few simple integrations, we arrive at: 



(Cc(xi,ti)Cc(x2,t 2 )) 



(19) 



i e 

648vr 2 e 



j5(n - t 2 ) sinc(ai/|xi - x 2 |) 



where sinc(:r) 
as 



_ sm(x) an( j gmc ^Q^ _ i us express C,' c 
C' c {x,T) = ^-G^r) X/ \{x,r), (20) 



36tt v 26' ' ^ ' n 
where the variable £ is a random Gaussian field with: 

< £(xi,ii)£(x 2 ,i 2 ) > = S(ti ~ r 2 ) sinc(ai?|xi - x 2 |). 

(21) 

Now, plugging the coarse grained value of £' from Eq. 
(|20|l into the left hand side of Eq. i|13|) . we ma Y obtain 
the coarse-grained field equation: 

,' = ^Me-^)*«M)]. (22, 



IV. THE CASE OF A QUADRATIC POTENTIAL 

In this section, we apply the coarse-graining method, 
developed in the previous section, to a chaotic inflation- 
ary scenario with a quadratic potential, and investigate 
the evolution of the background during inflation under 
the effect of second-order back reaction. 

Here, we assume that the potential for the inflaton 
field, (p, has a quadratic form: 



V(ip) = imV 



(23) 



Plugging this into Eq. (|10|) , and then lIl'L'l) , the equation 
of motion for reduces to 



9' + m 2 = — mQh(x,T). 

24-7T 



(24) 



This equation can be solved perturbatively, taking £ as 
the perturbation variable. can be approximated by 



= o + 0i + 2 , 



(25) 



where ©o, 0i and 2 are zeroth, first and second order in 
£, respectively. As the first step in solving Eq. pijl. we 
note that 0o is simply the expansion rate in the absence 
of perturbations, 



0o 



-m 2 r, 



(26) 



where our convention for r is so that it has a negative 
initial value Tj n and goes to zero towards the end of in- 
flation, so that decreases as inflation proceeds. Now 
we substitute ©o back into Eq. (|24|) to find 0i, 



(27) 
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This step is easily followed by calculation of 02- 

(28) 

Since £ is a random Gaussian variable, we are not in- 
terested in the spontaneous value of O, and instead in- 
vestigate the statistical behavior of (0), the expectation 
value of the expansion rate over time as more infrared 
modes cross the hubble radius and may build up a back 
reaction effect. As Oo does not have a £ term, and the 
expectation value of ©i vanishes over time, we have: 



(e> = e + (e 2 >, 



(29) 



where (82) can be derived from Eqs. I|28|) and l|21|l . which 
finally enables us to calculate expansion rate with correc- 
tions due to stochastic back reaction: 



(9) 



-m 2 T 



1728tt 



:(r 3 



J- 



(30) 



Contrary to our expectation from [5j, we can see that 
the back reaction term has a positive value and it is in- 
creasing with time, as if it is speeding up inflation. How- 
ever, as we see in the next section, this is merely a result 
of our choice to average over hyper-surfaces of constant 
proper time, while other choices of time hyper-surfaces 
may lead to different forms, or even a different sign, for 
the back reaction term. 

Nevertheless, as both and r are observable, we may 
conclude that the back reaction does not vanish during 
inflation and since its magnitude is increasing with time, 
it may become important if inflation lasts long enough. 
Notice that this correction is, in fact, due to wavelengths 
larger than the Hubble radius as, using Eq. (©2) 
could be expressed as an integral over k, the wave number 
instead of r, which would yield 



<© 2 ) = 



1 



594tt 2 



dink 



(31) 



and the above integral is only over modes which have 
already crossed the Hubble radius. In this sense, this is 
an infrared back reaction effect. 



Using Eq. J7J, we can see that ip may be used as a 
measure of e-folding: 



N e (r) =^i«-V(r), 



(33) 



and thus, in order to find how changes as a function of 
N e , it is sufficient to express it in terms of ip. Combining 
Eq. (jT§|l with £p — ip'-^p, we can re- write Eq. as 

where £ in the above equation like is, again, a random 
Gaussian field, defined by 

< ^1,^1)^2,^2) > = <5("0i - 1P2) sinc(ai7|x! - x 2 |). 

(35) 

In the case of single scalar field inflationary models 
with a power law potential 



(36) 



using the procedure described in the previous section, 
Eq. (|34|l can be solved to obtain in terms of ip: 



(6) =a a / 4 (3M 4 y/ 2 iP a / 4 + 

^ S - 3( ^ + a/ 2 _^ + a/a )<37) 



(3M 4 ) 3 / 2 a 1+3 «/ 4 



2887r 2 (a + 4) 

Note that we assumed ip ^ towards the end of in- 
flation. The result may be simplified for the case of a 
quadratic potential: 

(0) = V2(3Af 4 ) 17 V /2 

+ m ^f 2 r^-^ n ). (38) 

We see that, if averaged over the hyper-surfaces of con- 
stant e-folding rather than constant proper time, the sec- 
ond order back reaction has a negative sign, i.e. it slows 
down inflation. We again believe that this is a real effect, 
as both and N e are observable parameters. 



V. GENERALIZATION TO POWER LAW 
POTENTIALS 

In this section, we generalize the calculations of the 
last section to power law inflaton potentials. Also, rather 
than proper time, we use the inflationary e-folding num- 
ber as our time variable, as it has more direct and intu- 
itive observational relevance. The number of e-foldings 
since the beginning of inflation is 



N e (r) 



0dr. 



(32) 



VI. DOES BACK REACTION CHANGE 
INFLATIONARY PREDICTIONS? 

Let us calculate the ratio of the second order term in 
Eq. (|37|l to the zeroth order term 



(02 

0n 



3M 4 
288tt 2 



a 



(39) 



if tp -C ipin, close to the end of inflation. Currently, the 
most stringent constraints on the inflationary potential 
is due to the observations of the cosmic microwave back- 
ground anisotropies by the WMAP satellite 22]. The 
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FIG. 1: The example of an inflaton potential with a flat 
portion. Back reaction can dominate the slow-roll inflation in 
this portion of the potential. 



observed amplitude of anisotropics at the observable hori- 
zon scale require: 



4 tc 



k 3 P ( (k) 



2tt 2 



12tt 2 V' 2 ((p) 



= (2.1 ±0.3) x 1(T 9 , 



for k = 0.002 Mpc" 1 , 



with V' being the derivative of V with respect to ip. Then 
we can use the simple slow-roll prediction (first term in 
Eq. |2H1 along with Eq. I10|) to relate <p to xj), and express 
A 2 ^ in terms of ip: 



A/f 4 

n 12tt 2 y 



(41) 



Assuming that ipin — N e (Eq. I33fl . and that inflation 
ends when V ~ V (in Planck units), Eq. (|39|l can be 

expressed in terms of A-jz: 



<e 2 ) i 



-A^N e = 2 x 10~ 8 



2 x 10 



Although, this result is specific to the power law in- 
flaton potentials, we expect a similar order of magni- 
tude for any smooth potential that is consistent with 
observational bounds. Does this mean that back reac- 
tion will always be negligible during inflation? Of course, 
the number of c- foldings, N e , could be arbitrarily large. 
Indeed, after enough e-foldings, we end up with the self- 
similar chaotic inflationary picture which results from the 
stochastic formulation [23|]. However, only ~ 70 of these 
c-foldings could fit in our local Hubble patch today, and 
any effect on the power spectrum due to the back reaction 
of the super-Hubble modes is inaccessible to observations. 

Another possibility is that the inflaton potential may 
not be smooth. The observational bounds on A 2 ^ hold 
from the observable horizon (~ 10 4 Mpc) down to the 
scale of the Lyman-a forest (~ 1 Mpc) |24|, constraining 
only about 10 e-foldings in the inflaton potential. An un- 
usually flat portion in the inflaton potential (see Fig. 



may yield a significantly larger amplitude of fluctuations, 
e.g. A 2 ^ ~ 0.1, and be inaccessible to cosmological ob- 
servations if it happens at small enough scales. However, 
such a large amplitude may result in the production of 
Primordial Black Holes (PBHs), as the modes enter the 
Hubble radius, which results in further constraints on 
Ak, 

To study the production of PBHs [25j , let us look at 
the comoving scale of R < 100 kpc, with the primordial 
amplitude of fluctuations An. This scale enters the Hub- 
ble radius when aH ~ tt/R, which happens at redshift z, 
roughly given by 



10 9 i? _1 (kpc), 



(43) 



where Ho ~ 70 km/s/Mpc and fi TO ~ 0.3 are to- 
day's Hubble constant and density parameter respec- 
tively, while z eq ~ 3 x 10 3 is the redshift of radiation- 
matter equality. As the modes of scale R enter the Hub- 
ble radius, a fraction, /, of the Hubble patches have a 
fluctuation amplitude larger than 1, and thus may form 
PBHs, where / is given by 



(40) /(A R ) 



exp[ 



c 2 



2A 2 , J v /2^ 



R 



A TC exp[-A TC 2 /2]- 
(44) 

The current density in PBHs in units of the critical 
density, fipBH, is given by 



fi PBH ^/(A K )( — )e, 



(45) 



where 10~ 4 < e < 1, is the fraction of the mass in a 
Hubble patch that goes into the PBH [2(|. Requiring 
that fipBH < f^m yields 



f(A n ) < < i - 2 i?(kpc). 

ze 



(46) 



For A 2 n = 0.1, Eq. @£| gives f(A n ) 



10~ 3 , and thus 
Eq. 14t)|) . in combination with large scale structure ob- 
servations, require 



0.1 kpc < R < 100 kpc. 



(47) 



Other astrophysical observations may put tighter con- 
straints on a population of PBHs in the universe [27| . 
Instead of going through these constraints, let us give 
an example. Assume at R — 10 kpc, A^ = 0.1 and 
e = 10~ 4 . Using Eq. g5j), fi PB H ~ 0.01O m , while the 
PBH mass, A/pbh is eMjj, with Mjj, being the Hubble 
mass at the Hubble entry: 



H ~ 1O 17 M (^) 2 ~ (10 5 M )i? 2 (kpc). (48) 



M 



Therefore, with our choice of parameters, about 1% of the 
mass of the universe will end up in ~ 10 3 M Q black holes. 
This is consistent with current astrophysical constraints 
on such objects |27j . 
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To summarize this section, we showed that if we as- 
sume a smooth inflaton potential, in single matter field 
models infrared back reaction (during the last 70 e- 
foldings of inflation) should be negligible, at least to the 
2nd order (one loop) in the amplitude of perturbations. 
However, if we relax this assumption, it is possible to 
construct an inflaton potential that allows a significant 
back reaction effect during the slow-roll inflationary pe- 
riod, without breaking any observational constraints. 

VII. CONCLUSIONS 

We use a stochastic coarse-graining method to ap- 
proach the problem of infrared back reaction in a 
single scalar field inflationary cosmology. Using the 
coarse-grained generalized Bardeen parameter, we self- 
consistently derive the stochastic Langevin equation that 
governs the evolution of the local expansion rate. Solv- 
ing this equation to second order (one loop) for chaotic 



inflation with a power law potential, we see that there 
is a non- vanishing average back reaction term due to the 
infrared modes that have left the Hubble radius. Al- 
though, the size of this effect will be negligible (within 
the observable universe) for smooth inflaton potentials, 
we show that it is possible to have non-smooth poten- 
tials which allow significant back reaction happening on 
sub-horizon scales, without breaking any observational 
constraints. A possible implication of such non-smooth 
potentials is the generation of a population of massive 
primordial black holes, which could be further used to 
constrain or verify these models. 
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scalar field in a de Sitter back ground. 
The cosmological scale factor is denoted by a(t), the Hub- 
ble expansion rate by H , and the scalar matter field by 
<P- 

For a more detailed discussion about calculation of this 
parameter see |gL 1 1 1| . 

In linear perturbation theory, 5(j> is forced to be equal to 
Sip by the off-diagonal dj equations for all theories in 
which the off-diagonal components of Tij vanish to the 
linear order. This is in particular the case for scalar field 
matter. 

This is the simplest example of a cutoff function. In gen- 
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eral, one can use any function which is close to one for to eliminate the ultraviolet divergence. 

k < aH and tends to zero as k goes to infinity fast enough 



